SCHWINGER'S FORMULA AND THE PARTITION FUNCTION 
FOR THE BOSONIC AND FERMIONIC HARMONIC OSCILLATOR 



^1- 

On 



L. C. Albuquerque, C. Farina * and S. J. Rabello 



5_i . 

Instituto de Fi'sica 
Universidade Federal do Rio de Janeiro 
(N : Rio de Janeiro, R.J. CEP - 21.945-970 - Brazil 

> 
(N 
O 
O 

m 
o 
^t- 



a: 

^ . Abstract 

> : 

•rH , 

We use Schwinger's formula, introduced by himself in the early fifties to compute 
effective actions for QED, and recently applied to the Casimir effect, to obtain the 
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The computation of the partition function of the usual harmonic oscillator is 
probably one of the most elementary exercises in a Statistical Mechanics course. 
There are many ways of making this calculation, and undoubtedly the easiest is the 
direct one, that is, 

Z{(5) :=Tre- pft 



= E 



-/3(n+±)u, 



e ^""2, 

n=0 

e 2 
~~ 1 - e-P" 

= 1 ~1T- . (!) 

2sinh(^) 

where we simply used that the eigenvalues of the Hamiltonian operator for the har- 
monic oscillator are given by (n + \)uj, with n = 0, 1, 2, ... (we are using % = 1) and 
summed the infinite terms of a geometric series. 

However, it is exactly the simplicity of handling with this example that makes 
it a perfect "laboratory" to test or develop other methods of computation, as for 
instance, the path integral method^ 1 -*, the Green function method 1 - 2 ), etc.. 

Our purpose in this letter is to apply a formula invented by Schwinger in 1951 ( 3 ) 
to compute effective actions for QED, and recently applied with success in the com- 
putation of the Casimir energy for both the massless and massive scalar field^ 4-6 ), 
to obtain not only the result of equation (1), but also the partition function for 
a fermionic harmonic oscillator. Curious as it may seem, this approach has never 
appeared in the literature. 

Let us start with the bosonic case. It is well known that the corresponding 
partition function can be written as^ 7 ) 

Z(/3) = det -Hu 2 -d 2 T )\ Tp , (2) 

where the subscript T v means that the operator u 2 — d 2 acts only on a set of functions 
which are periodic, with period (3. In ref. [7], Gibbons used the generalized £- 
function method to compute such a determinant. Here, we shall use Schwinger's 



formula (deduced in the Appendix A): 



1 



\nZ((3) = -Tr / dss 



,y—\ —isL u 



(3) 



where L u := u 2 — d 2 . and we chose the regularization based on the analytical contin- 
uation, instead of Schwinger's one. 

For periodic boundary conditions we have 



Tre 



-isL u 



= E 



-is[" 2 +n 2 (%) 2 ] 



(4) 



so that 



-I poo 

1 n= _oo^0 
1 oo 

= l^F(u) + 

where we introduced the one-dimensional inhomogeneous Epstein function 



2tt 



T{y)E» (u,l), 



(5) 



(6) 



n=l 



defined [i = |^ and used the well known integral representation of the Euler Gamma 
function J °° dtt"- 1 e~ at = a~ v Y{v). 

Although the above series converges only for Re v > | , it can be analytically 
continued to a meromorphic function in the whole complex plane given by( 8 -* 



E? (u, 1) 



i + v /^i>-±) 



2^fu ^ f n7r 



K u _i (2vrn/i). 



(7) 



It is worth noting that the structure of poles of E± (V, 1) is governed by the 
poles of T{v — |). Hence, they are located at v = |, — |, — |, and so on. As we 
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see, E± (u, 1) is analytic at the origin. 



Substituting (7) into (5), taking the limit v — > and observing that the divergent 
terms cancel without any further subtraction, we get 



n=l 



(8) 



In order to compute the summation on the r.h.s. of (8), we appeal to the 
formula^ 9 -* 

Ki {nu(3) = 



(9) 



2nu(3 

Inserting (9) into (8), and using equation (B.l) (see Appendix B), we obtain the 
final result 



InZ(P) = - In 



2 sinh 



u/3 



i 



2 sinh 



(10) 



in perfect agreement with (1). 

For the fermionic case, it can be shown that 

Z'(f})=det +1 (u 2 -d 2 T )\r a , 



(11) 



where the subscript jF a now means that the operator acts only on a set of an- 
tiperiodic functions. As a consequence, the eigenvalues turn to be X n = u 2 + (p 1 ^ , 
with p an odd integer. Then, using these eigenvalues in the analogue of equation (3), 
but remembering that for the fermionic case, instead of the factor | we must write 
(-1), we get 



\nZf(P) = -T(v) 

p=odd 



7T 



Adding and subtracting T(u) J2 p =e 
tion, we may write 



(12) 

to the r.h.s. of last equa- 



lnZ f (P) =2T(u) 



2v 



(13) 
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where now ji = Following exactly the same steps as before, that is, using the 
analytical continuation of the one-dimensional inhomogeneous Epstein function given 
by (7) and taking the limit v — > 0, we get 

lnZ'(/?) =(3u + 4^ f; f-K_ h {nuP) - 4^ £ {^K.tfw®. (14) 

n=l n=l 

Using equations (9) and (B.l) we finally obtain, 



\nZ f (/3) = In 



-l 2 

sinh (uj/3) 



sinh ( ^ 



Z f (/3) =4 cosh 2 (^pj . (15) 



This result can be checked easily in the following way: the fermionic oscillator 
that is being considered here is the second order Grassmann oscillator studied by 
Finkestein and Villasante* 10 ) (in fact, we are dealing here with the particular case of 
iV = 2 of their work). In this case, it can be shown that there are only three energies: 
(double degenerated), 
+(3uj, —flu, so that, if we use directly the definition for ((3), we will obtain 

Zf(0)=Tre-P AI 

= 2 + e? u + e"^ 

= 4cosh 2 (^). (16) 

In this paper we applied the Schwinger's formula for the one-loop effective action 
to the computation of the partition function for both the bosonic and Grassmann 
harmonic oscillator. Regularization by analytical continuation was adopted. We note 
that depending on the boundary condition choice ( eg. Dirichlet condition), new 
subtractions (renormalizations) can be needed. But this is easily done remembering 
that the Schwinger's formula contains an integration constant, which can be used to 
subtract these divergent terms. 

One of us (CF) would like to thank M. Asorey, A. J. Segm-Santonja and M. 
V. Cougo Pinto for helpful discussions on this subject. This work was partially 
supported by CAPES and CNPq (Brazilean councils of research) . 



Appendix A 

Let 

e T ^ = detL^ = exp Tr In L u \r, (A.l) 

where L w = uj 2 — d 2 , and the subscript T means that some boundary condition is 
assumed. Hence 

I» = Tr \n(u 2 - d 2 T )\ T (A.2) 

Taking the variation with respect to uj 2 (the subscript T will be omitted but the 
boundary condition is understood), we obtain 



(A3) 



POO 



dse~^ L --^ 5L U 



= S 



DC 



ds 



Tr / — e - i ( L ^- i < E ) 



o 



Apart from an additive constant to be fixed by normalization, integration leads to 

ds 



-Tr [ 
Jo 



T(u) = -Tr I —e- tsL ", (AA) 



where the ie factor of convergence is tacitly assumed. Last equation is clearly ill- 
defined. One way to circunvent this problem was pointed out by Schwinger: one 
introduces a cut-off s , so that the integral can be performed. Then, subtracting the 
divergent terms for sq = 0, the remaining integral is finite. However, there is another 
possible choice^ 5 ), which consists in replacing expression (A. 4) by 



\u) = -Tr [ 
Jo 



oo 

v — \ —isL 



dss v - l e- lsL % (A5) 



with v big enough, such that integral (A. 5) is well defined. Then, after the integral 
is made, we perform an analytical continuation to the whole complex plane. After 
subtracting the poles at v = (when they exist), we take the limit v — > 0, and this 
way we get a finite prescription for the integral (A. 4). In the main text we adopt this 
approach instead Schwinger's one. 
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Appendix B 

In this Appendix we shall prove that 



-e~ na = - - In 
n 



2 sinh 



With this purpose, we will define S(a) such that 

S(a) := J2 -e~ 



a 
2 Jl 



a > 0. 



(B.1) 



(B.2) 



n=l 



Diferentiating both sides of (B.2) with respect to a and using the well known result 
of the sum of the infinite terms of a geometrical series, we get 

dS(a) 



e a 



da ' 2 sinh (f) ' 
Now, integrating in a we readly get 



(5.3) 



a 



a 

2 sinh — 
2 



(5.4) 



where we used that S(oo) =0. This completes the desired proof. 
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